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Abstract: We present here results on the fine structure of the static qq potential in d = 4 
SU(3) Yang-Mills theory. The potential is obtained from Polyakov loop correlators having 
separations between 0.3 and 1.2 fermi. Measurements were carried out on lattices of spatial 
extents of about 4 and 5.4 fermi. The temporal extent was 5.4 fermi in both cases. The 
results are analyzed in terms of the force between a qq pair as well as in terms of a scaled 
second derivative of the potential. The data is accurate enough to distinguish between 
different effective string models and it seems to favour the expression for ground state 
energy of a Nambu-Goto string. 
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1. Introduction 

The string picture of hadrons jl], ||, has a long history starting from the 1960's. Although 
this picture could explain the Regge trajectories of hadrons, it required 26 space-time 
dimensions to be a consistent theory Q. Effective string theories which could be quantized 
in any dimensions were introduced by Polchinski and Strominger in ||. In QCD such 
strings have the interpretation of gluonic flux tubes between a quark and an antiquark. 
Nambu made formal connections between QCD and the string models || . 

A long distance 1/r term in the quark-antiquark potential, distinct from the short 
distance couloumbic 1/r term was first observed in 0. Subsequently Liischer showed the 
universality of this term and since then it is known as the Liischer term ||, having the 
value —7r(d — 2)/24r where d is the number of space time dimensions and r the length of 
the string. 

In the eighties || and |i0| ] claimed to find the Liischer term in d = 3 SU(2) and 
d = 3 + 1 SU (3) lattice gauge theories respectively. However precise data on large Wilson 
loops or Polyakov loop correlators, necessary to identify the class of the effective string did 
not exist at that time. 

In recent times, significant progress has been made due to increase in computing power 



as well as improvement in algorithms. See flT| , 12, 13, 14, 15, 16] for example. 

Related to the effective string descriptions is the issue of the spectrum of string exci- 
tations. Recently analytical studies, based only on symmetry principles have been carried 
out in |17], [18|| . Extensive numerical studies have been carried out in [19] and 13, 14]. We 



do not discuss this topic here. One can refer to the review article by Kuti ]^0| for a general 
introduction and more details. 

In this article we present results of our simulations of the Polyakov loop correlators 
for d = 4 SU(3) Yang-Mills theory and compare the resulting static potential with both 
perturbation theory and string model predictions. A preliminary report of this simulation 
was presented in [21]. 
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2. c(r) at short and long distances 



In this article we look at c(r) = ^r 3 F'(r) where F(r) is the force between a quark and 
an antiquark. This quantity asymptotically tends to the Liischer term for large r. We 
measure c(r) at intermediate distances from slightly below 0.5 fermi to about 1.2 fermi on 
a 32 4 lattice at (3 = 5.7. It is therefore interesting to compare the values one obtains from 
perturbation theory at short distances as well as from the string picture. 



c(r) can be easily computed from 2-loop perturbation theory using the results in [22 
and p3fl. We just quote the final result for the sake of completeness: 



r 



c( r ) = -C F {aMs + c o( a Ms) } _ o d r a MS I 1 + 2c o«ms) 



+ .... (2.1) 



Here Cp = | and cq = 87r/?o (7 — §§)• ^he function a^jg is given by [23] 



with fi = — ln(r 2 A|jg) and fi = ln(/i). For SU(3) the constants /?o, 0i and 7 have 



their usual continuum values of , and 0.57722 respectively. The result of this 

computation is plotted as the curve '2-loop perturbation theory' in fig with the upper 



and lower curves corresponding to estimated upper and lower limits of Aj^g[22] 



At the other end of the length scale one expects a string like behaviour. Since non- 



bosonic forms of the string have been ruled out [24], our case of interest is the Nambu-Goto 



string. The potential in this case was first given by Arvis |2q] to be 



VArvis(r) = \/a^ - (2.3) 

An interesting point, noted even before Arvis is that the qq potential becomes purely 
imaginary for r < r c where r c = \ |26||, which is around 0.3 fermi. This behaviour 



was connected to the tachyon instability of the Nambu-Goto string [27|. It is curious that 
perturbative calculations also break down at about the same distance. 

Another form of the potential that is of interest is the so called truncated Arvis po- 
tential which is obtained by expanding the Arvis potential in a power series and retaining 
the first three terms: 

(rf-2)7T (d-2) 2 7T 2 

3. Simulation 
3.1 Algorithm 

Since we are interested in ground state properties, we measure Polyakov loop correlators 
as it is well known that the Polyakov loop correlator for large temporal extent strongly 
projects onto the qq ground state. We use the multilevel technique of Liischer and Weisz 
for measuring the Polyakov loop correlators very accurately. For details of the algorithm 



we refer the reader to [12] 
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This algorithm has several optimization parameters and the most important among 
them seems to be the number of "measurements" used to compute certain intermediate 
expectation values. We will refer to this number henceforth as "iupd" . Another parameter 
is the thickness of the time slice. For our case a thickness of two was optimal. 

It is also well known that improved observables can be constructed by replacing a bare 
link by its group average keeping the environment unchanged. This is known as multihit 
and we employ that too on the time- like links for the correlators. The group average is 
computed through the single-link integral defined by 

(U) = Z(J,J^r l dZ i J ;f ] where Z(J,J*)=[ [dU]e^ u Jt +ut J \ (3.1) 

OJ J JSU(3) 

For SU (3) the group average cannot be easily carried out analytically and is most often 
evaluated using a monte-carlo method. An alternative semi-analytic method for this aver- 



aging was proposed by de Forcrand and Roiesnel [28]. We used the semi-analytic method 
for the multihit and we estimate that this resulted in a 60% speedup of the code compared 
to using the monte-carlo method to reach similar levels of accuracy. 

In our test runs, for Monte-Carlo multihit the optimal value for iupd occurred at the 
same place for the two different values of r we considered but not for the semi-analytic case. 
To obtain maximum gain from the multilevel scheme, we tried to choose iupd close to the 
optimal values for our largest r. Another interesting observation is that the penalty paid 
for operating at non-optimal values of iupd seems to be higher for Monte-Carlo multihit 
than the semi-analytic method. 

We have carried out simulations at (5 = 5.7 on both 24 3 x 32 and 32 4 lattices using the 
Wilson gauge action. The lattice spacing at this f3 is 0.17 fm so that the temporal extent of 
the lattice is 5.4 fm while the spatial box is (4fm) 3 in one case and (5.4 fm) 3 in the other. 
Details about setting the scale is discussed in the next sub-section. 

For separations r = 2 — 6 each measurement involved simulations on 24 3 x 32 lattices 
with iupd= 12000 and 500 measurements were made in all. For the larger separations 
r = 5 — 9 simulations were done on 32 4 lattices with iupd=48000 and in all again about 
500 measurements were obtained. 

We also found that for the larger separations (r = 6 and higher) , we had to go to larger 
lattices to continue to gain from the multilevel scheme. On the 24 3 x 32 lattice increasing 
"iupd" even by an order of magnitude did not seem to help reduce the error on the r = 6 
correlator. However going to a larger lattice (32 4 ) did help significantly. 

3.2 Results 

In our simulations we measure the Polyakov loop correlator {P*P)(r) where r is the sepa- 
ration between the two Polyakov loops. The qq potential is determined as 

V(r) = ~log(P*P)(r) (3.2) 



where T is the temporal extent of the lattice. In our analysis we look at the force between 

IV (r 
dr 



the gg-pair given by F(r) = dV }^ and the scaled second derivative c(r) given by c(r) 
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Table 1: Results: Polyakov loop, Force, c(r) 



IT d • ^' ie ^ a ^ter tends to the asymptotic value of — for large r which is nothing 

but the Liischer term. We are interested in how this quantity approaches its asymptotic 
value. On the lattice these quantities are defined by 

F(f) = V{r) - V(r - 1) (3.3) 
c (f) = T —[V{r + 1) + V(r - 1) - 2V(r)) (3.4) 

where f and r are defined as in |^] to reduce lattice artifacts. Our results are presented 
m table [TJ. 

Since different r values are measured in the same simulation, they are quite strongly 
correlated. Thus when one evaluates the differences in potentials, it helps if one evaluates 
the differences for each measurement and then averages over different measurements rather 
than the other way round. In practice we compute the difference between jackknife bins and 
then compute the jackknife error for the differences. We also checked that our combination 
of one heat-bath with three over-relaxation sweeps and restricting the measurement to 
every fifth update resulted in negligible autocorrelation. 

The qq potential contains an unphysical constant which masks the properties of the 
flux tube to a certain degree. We therefore look directly at the force defined as in eq. (|3~3|). 
In figure [l] we have plotted the force versus 1 /r 2 . We use the force to set the scale on the 
lattice via r^Ffa) = 1.65 where ro is the Sommer parameter which we take to represent 
the physical distance of 0.5 fm. In our case the Sommer parameter turns out to be 2.93 
lattice spacings. 
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Table 2: Fit results 
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Figure 1: Force in d — 4 SU(3) case. The inset shows the last three points in greater detail. The 
fits 1, 2 and 3 are to a straight line, truncated Arvis and full Arvis respectively. Details about the 
fits are given in table 0. 



We do three different fits to the force data whose results are shown in table |2| . First 
we fit to the form sx + c where x = 1/r 2 with s and c as fit parameters. From the intercept 
of this fit we obtain the string tension to be aa 2 = 0.1554(1) . The slope of this line gives 
us the effective c(r) over this range of r and this turns out to be 0.304(3) which is still 
16% away from the asymptotic value of 7r/12. This clearly shows that we are still some 
distance away from the region where the model independent leading order behaviour is all 
that matters. 

Next we fit to the form s + cx + 3c 2 x 2 /2s which is inspired by the truncated Arvis 
potential eq. (|2,4|) . This fit gives us the string tension via s as 0.1559(1) and the effective 
c(r) to be 0.266(3) which within errors is almost identical to the universal value of 7r/12. 
Finally we fit to the form expected from the full Arvis potential. Within the fit range and 
errors, this fit is indistinguishable from the fit to the truncated Arvis potential. Unfortu- 
nately we cannot yet do an analysis of the the type envisaged in [17] where an additional 
coefficient in front of the 1/r 4 term can be determined by fits. Such an analysis would 
probably require data with significantly reduced error bars. That seems to be out of scope 
at the moment. What we do see is that the 1/r 3 term in the potential for open strings 
with fixed end boundary conditions is indeed consistent with the Arvis potential and any 
correction to this would have to be really small. 

In figure ^ we have plotted the scaled double derivative c(f) as a function of r/rQ. The 
horizontal line at c(r) = — 7r/12 is the asymptotic value. Also shown are the predictions of 
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Figure 2: Scaled second derivative for d = 4 SU (3) case. The circles are from p9| 

c(r) from the Arvis as well as the truncated Arvis potentials. The perturbative predictions 
are plotted for small values of r extending upto about 0.25 fermi. 

4. Discussion and conclusions 

In this article we have looked at how the QCD-string behaves at intermediate distances. The 
only other simulations of the Polyakov loop correlators in d = 4 SU (3) theory are reported 
in ^9|. To compare our results with them we have plotted their data along with ours 
(circles in figure ^). In the range where we overlap, the agreement of the data is excellent. 
We have been able to extend the range modestly albeit at a considerable computational 
cost 1 . However even this extension is important as our two extra points definitely show 
that the data falls quite nicely on the Arvis curve which was not so clear yet with the 
previous simulations. We unfortunately cannot really distinguish between the truncated 
Arvis potential and the full Arvis potential as the difference between the two curves is less 
than 1% already at 2ro or 1 fermi. 

Polchinski and Strominger || found the Liischer term as the leading 1/r correction 
in their effective string theory. It has now been claimed (for closed strings) p0| , [H], 32] 



that even the coefficient of the 1/r 3 term is the same as the one obtained from the Arvis 
potential. We find it interesting that the same seems to be true for the open strings also. 
While in effective string theories there is no apriori reason to expect the full Nambu-Goto 
behaviour, at the moment we do not see any deviations from it beyond a certain distance. 



1 In it was estimated that for a given lattice, the computational cost to reach the same level of 
accuracy in c(r) grows as r 4 with r. 



-6- 



This simulation was performed on a comparatively coarse lattice. Existing results on 
finer lattices [29] have shown, for the smaller r values, that the data lies above the values 
obtained in this simulation. If that behaviour persists for larger r values too, then the scale 
at which convergence with truncated Arvis potential sets in would be larger. It would be 
really interesting to study the continuum limit of this scale. However that is work for the 
future. 
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